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FOREWORD 


This report explores the effect of fog and cloud droplets on aircraft wake 
vortex decay rate. The question was originally posed by Dr. Vern Rossow of 
the NASA Ames Research Center. The report is exploratory and should be 
considered preliminary. 

The research was conducted by Dr. T. H. Moulden and Dr. Walter Frost 
at the University of Tennessee Space Institute under the technical direction of 
Dr. George H. Fiehtl of the NASA Marshall Space Flight Center. Resources 
for the conduct of the effort were made available by Mr. John Enders of the 
Aeronautical Operating Systems Division, Office of Advanced Research and 
Technology, NASA Headquarters. 
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I. 


Introduction 


It has been suggested that the trailing vortex system 
behind an aircraft may decay more rapidly when fog is present 
in the atmosphere. Such an effect is of significance to 
airport scheduling where the hazards to aircraft operation 
associated, with the strong vortices generated, by the current 
class of large jet aircraft require time consuming separation 
distances between successive landings and. take-offs. As 
airports become more congested., this situation creates flight 
delays and financial penalties. Thus, there is ample need 
to study all possible conditions under which the rate of 
vortex decay is increased.. 

The mechanism by which fog droplets would, alter the 
rate of vortex decay is not fully understood.. One can envision 
that the acceleration of the water droplets as they first 
enter the vortex would, remove energy from the flow and. retard, 
its circumferential velocity. Once the droplets are 
entrained, into the vortex environment they will absorb very 
little additional energy — save for a redistribution of the 
vortex velocity profile as particles near the center of the 
vortex migrate towards the outer reaches. This particle 
motion will be discussed in detail later. Due to the intro- 
duction of the water droplets, the fluid has assumed a two- 
phase character and. may well suffer some change in apparent 
physical properties' 1 ^. Thus, an equivalent fluid viscosity, 
taking into account the influence of particles on the turbu- 
lent eddy viscosity, could, be defined which, intuitively 


suggests a retardation of the vortex, but physically may 
produce higher Reynolds numbers thus giving greater inertia 
to the two-phase flow and thus sustain the vortex life. 

Another effect that may operate the phenomenon of 
(2 3 ) 

vortex bursting * . Due to the changes in physical proper- 

ties and in the velocity profiles brought about by the water 
content, it is possible that vortex bursting would occur 
prematurely. As the understanding of vortex break-down 
is not well documented, very little can be said on this 
point at present, but stability theory does point out the 
significance of the local velocity gradients. It is, however, 
felt that a solution of the full Stokes-Navier equations will 
be necessary before this problem of vortex bursting can be 

(4 ) 

fully appreciated. . It is urged, that such a study be 
undertaken. 

In the present report, some comments are made concerning 
the motion of non-interacting spherical particles in a 
vortex-type flow. The remarks are made in the context of 
an exploratory exposition rather than of a definitive study. 
Engough is said., however, to clarify the direction that future 
studies should, take. 
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II. Preliminary Remarks 


Before embarking upon the study in hand, a remark 

concerning the nature of fog is in order. 

( 5) 

According to Tverkoi v ' the distribution of water 

droplet radius in a fog is characterized as in Figure 1. 

-4 

Typically, the mean particle has a radius a=8 x 10 cm and 

-9 

a mass m = 2 x 10 gm. It will be assumed in the present 

study that this median value is independent of the liquid 

water content present in the atmosphere. This may not be the 

case when the atmosphere is very near to its dew point. 

When free to fall in a still atmosphere, this particle 

of mean diameter will have a motion described by the equation: 

du 6Trauu (1) 

dt -g “ m 

In this equation, u is the velocity of the particle at time 

t. The fluid viscosity coefficient is denoted by (x. It 
is assumed in writing equation (1) that the particle exper- 
iences a retarding force given by Stokes* formula v ' for 
the slow motion of spherical bodies. Due to the smallness of 
a, surface tension dictates that the particle is essentially 
spherical. Equation (1) readily solves to show: 

- 6Trixat jj (2) 

for a particle at rest at t = o. This velocity is plotted 
as a function of time on Figure 2. Also included on this figure 
(for comparison with later studies) is the numerical solution 
to equation (1) obtained by means of a second order Eunge-Kutta 
technique. 


u = 1 1 -exp 
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The terminal velocity for this mean diameter particle 
is 0.81 cm/sec giving a terminal Reynolds number of 0.01. 

The Stokes flow assumption is thus seen to be in order. 

The adequacy of the numerical solution is also evident from 
the plot of Figure 2 and. this observation is pertinent to 
the subsequent work. 

If the characteristic velocity u c = mg/6irna and characteristic 

2 

time t c = m/6ir^ia are defined then the Froude number F^ □ u c /ga 
and the Reynolds number Re=u ap /\x relate by the equation 
Re = 9pF,/2p . Equation (2) indicates that when t = t then 

Q. W C 

2 

u is approximately 33% of u c » Clearly, t * F^a /Re 7 so 

2 

that for fixed physical conditions in the problem t ~a . 

c 
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III. Forces on a Spherical Particle 


Figure 3 illustrates the coordinate system employed and 
the forces acting on the particle. Four forces will govern 
the motion of the particle. These include the gravitational 
force, the lift, drag and centrifugal forces. From the onset it 
is assumed, that the fluid motion is laminar in nature so that the 
influence of a turbulent structure on the particle motion need 
not be considered.. It is further assumed that the particle 
is spherical and of uniform composition. Because of this 
assumption, the particle motion may be treated as that of a 
point mass. The particle will experience no torque. Despite 
the fact that the local velocities inside the vortex may be 
considerable, it will still be assumed, that the drag force 
experienced, by the particle is given by Stokes' formula. 

As this equation is lineal* in the relative velocity between 
the sphere and. the fluid, we are at liberty to resolve the 
drag force simply by using the velocity components. Possibly 
during the initial phase motion, when the relative velocity 
between fluid, and particle is large, a more accurate expression 
for the drag force would, be in order. However, in view of the 
large simplifications involved, with the assumption of laminar 
motion, it is felt that such an elaboration would be misleading 
and inappropriate. 

In this study it is assumed that the vortex axis is hori- 
zontal so that the gravitational force which acts downwards is 
perpendicular to the axis of the vortex. 
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Due to the curvature of the streamlines within the vortex 
there is a force induced on the particle which is directed 

( 7 ) 

towards the center of the vortex. Within first order accuracy' ' 

3 2 

this force is stated as 2irpa w /r. Here w is the relative 
velocity component in the asimuthal direction and r the 
distance of the particle from the vortex center. Finally, 
there is a centrifugal force acting radially outwards. 

It is useful to compare the magnitudes of the various 
force components. To do this we note that the force components 
can be written: 

F_ = X,(U-u) 

z ■*- 2 

F r = V V " u) ~T"'< w ~ w)2+ g 1°S 9 + F (3) 

F = X (W-w) -g 
6 


where 


_ 631 a _ 3 Q_ 

K 1 ~ m » k 2 ~ 2 pw 


For the standard size particle it follows that: 

X x = 1.2. x 10 3 sec - 1 ;X 2 = 0.00559 


Thus with the exception of positions very close to the vortex 
core (r«l) the drag force predominates the motion. As a 
result it is expected that the particles will accelerate rapidly 
to the local flow condition and then move essentially with the 
streamlines of the flow. 

It should, also be pointed out that the particle mass 
only enters equations (3) via the drag term (X^). Here 
the relationship is an inverse one so that the heavier particle 
responds more slowly to the fluid, motion. 


6 



Of course in the real flow where the motion is turbulent, 


these conclusions would, need some modification. 
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IV. On a Model for the Vortex Flow 


Let w =(£, T |, C ) be the vorticity vector in cylindrical 

coordinates (z,r,9). If in addition, the velocity components 

are (U,V,W,) then for a flow that is axially symmetric, given 

by the condition d/d© =0, it follows that : 

« 1 dWr 

q r dr 


_-l dWr 
^ r dz 


. dV dU 

c= S ' ~5T 

The continuity requirement 

d rU j_d rV = o 
+ ~dF~ 


is satisfied by a streamfunction Y such that 


U 


1 dY 

r dr 


V = 


-1 dY 
r dz 


The present model of a line vortex will admit variations 
in all three velocity components (U,V,W) but makes the 
assumption that the vorticity vanishes in the (r,z) plane. 
This demands that £= 0 ; then we find : 


=0 


(4) 


d 2 X , 1 dX , d 2 X X 

5* + r 3F r * 

where X = f/r denotes a modified streamfunction. 

If it were further assumed that both £ and q identically 

vanish, then it would, follow that Wr= const; in conformity 

with the potential vortex flow. However, to remove the 

(9) 

singularity on the axis, the result given by Lamb 


Wr 


-z?{ 1 - exp 


-r 2 Uoo l] 
4vz J / 


(5) 
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has been adopted.. Away from the vortex core, the flow is 
again essentially irrotational. 

In this model of the vortex flow it is interesting to 
note that the flow in the (r,z) plane (i.e., the U,V 
velocity components) is completely uncoupled from the motion 
in the asimuthal direction. There is, then, no mechanism 
for the singularity on the axis, that is evident from the 
solution w <v 1/r, to be removed by an axial velocity component 
along the axis^°\ This latter coupling is the province 
of the viscous nature of the flow: which is ignored, herein. 

It was suggested by Moore and Saf fman^^and confirmed 
experimentally v , at least for one value of Reynolds number, 
that the axial velocity along the core of a trailing laminar 
vortex can be expressed as : 

= ^ g - { 1 + 8 . 57 o 2 */5i~x 10~ 5 ) (6) 

y/z/c 

where a is the wing incidence in degrees and c the wing chord 

length. Re is the Reynolds number of the flow based, upon 

the wing chord and. is sufficiently small for the vortex flow 

to be laminar. Hence, far from the wing it is expected, that 

the axial velocity will decay as the inverse square root 

( 12 ) 

of the axial distance. The data shows that equation (6) 
loses its validity 'close' to the wing (z/c<100). 

In the present study, equation (4) is solved under an 
assumed initial distribution of axial velocity. The initial 
value to be adopted, for u/U is still in question but could 
be selected so as to satisfy equation (6) in the far field.. 
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Using a Fourier-Bessel expansion, it is readily found that 


U(r,z) = 2 a J (ix r) exp (-^ z) 
n=l ° n n 


(7) 


V(r,z)= 2 A n r) exp (-M^z) 

n=l 


¥(r,z) = r 2 An J, (|a r) exp ( — j_i z ) 

__ i M *■ fl U 

n=i ^n 


where 


3>(r,z)= -2 An J, (|x r) exp (-a z) 
n=l (xn A n n 


A n 


— ~ k- \ r U (r) J (m- r) 

RfJ^^R )] 2 J ° 


dr 


with ll the zeros of J (u R) = 0 
n o n 

and U Q (r) is the assumed initial axial velocity distribution. 


An expression like 


U Q (r) = 


l+k 2 r 2 


( 8 ) 


has the correct form for specified arbitrary constants k and k^. 

The far field boundary conditions have been satisfied on 
r = R where R can be several wing chord lengths in magnitude. 
Further comments concerning the solution of equation (4) are 
contained, in the Appendix. 

It is apparent from equation (7) that the axial velocity 
component along the axis dies away exponentially far from the 
origin of the vortex. This results contrasts with that 
contained in equation (6). A typical calculated distribution 
of axial velocity component along the axis is shown on Figure 4. 
Also plotted, on this figure is the variation of axial velocity 
given by the approximation presented, in the Appendix (equation A3) . 
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The differences between the exponential and. inverse square 
root decay rates is an indication of the fashion in which the 
solution is influenced by the inclusion of a viscous term. 

Typical calculated axial and radial velocity components 
are shown in Figures 5 and. 6 respectively. 
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V. 


An Indication of Particle Motion 


As a prelude to embarking upon a discussion of the 
particle motion within the vortex, some orders of magnitude 
of the relevant terms can be written down. 

First, a comparison can be made between the response of 
the particle to a change in the fluid velocity and the 
influence of gravity. For a particle subject only to a 
Stokes drag force, the equation of motion states 

it - sj* <o.-u) <9) 


if is the velocity of the fluid, environment, then easily 



1 -exp 


- 6nnat 

m 


which contains an identical exponential term to that in 
equation (2). Hence the response to gravity and changes in 
fluid velocity are enacted in the same time scale. This, of 
course, is clear from the physical nature of the situation 
as embodied in the respective differential equations. 

It was seen from Figure 2 that the time required for 

_3 

u/U ~ 0.99 is t = 4 . 6m/6im.a seconds or t = 3.82 x 10 sec 

00 

for the standard particle. This suggests, as was noted, above, 
that the general mode of motion will be for the particle, 
initially at rest, to very rapidly accelerate to the local 
fluid velocity and then progress with the fluid environment. 
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Hence, work done on the particle will predominently be over 
the very short initial time period associated with the rapid, 
acceleration phase of the motion. After this, the particle 
motion will deviate little from that of the fluid, and. the 
rate of work is extremely small. 

To describe the motion in more detail, the equations of 
motion (3) were solved, using a second order accurate Runge- 
Kutta integration procedure. The accuracy of this solution 
algorithm as applied to the current type of problem is shown 
in Figure 2 and. was discussed previously. The velocity field 
U, V, W of the vortex was taken from the solution discussed 
above. Then with F , F and F~ as in equation (3) we have 

Z 1 v 


du 

dt 


dv 

dt 



dw 

dt 



( 10 ) 


u= 


dz ; 
dt 



w = r 


d£ 

dt 


which give the particle position. 

The computer program developed, to solve these equations 
returns the particle position, velocity components, accelerations, 
and the total work done on the particle as a function of time. 
Figure 7 shows typical traces of the work done on the particle. 
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As expected from the tentative discussion of the nature of the 
motion, the rate of work is very large over the initial time 
period while the particle attains the fluid velocity. After 
this, very little work is performed. The traces of Figure 7 
are for particles of the same mass and. starting from different 
asimuthal positions at the same radial distance from the 
vortex axis. It is seen that, to first order, there is little 
influence of asimuthal starting position on the work done. 

Figure 8 shows how the motion of particles of various 
masses progresses with time. All particles start at the same 
point in the vortex. 
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VI. An Estimate of the Total Energy 


Let 6W be the work done on a single particle. If there 
are N particles in a volume dV then the total work is 

W = f N SWrdrdOdz 
dv B t 

=2irNU / f r Swdrdt 

oo J o J O 

as pointed out before, the integral in time need, only be 
carried out over the initial acceleration phase of the motion. 

In evaluating this work integral it has been assumed, 
that the water content of the air is contained as a uniform 
distribution of mean sized water droplets. While not realistic, 
in light of the data presented in Figures 1 through 8, it 
is adequate for an order-of -magnitude calculation. It is 
also assumed, that the work done is independent of asimuthal 
position (see Figure 7) so that one integration can be carried, 
out analytically. 

(5) 

Using a typical value of liquid water content' ' of 

3 

0.2 gm/m , evaluation of the above integral gives 

W/U = 0.0022 Newton sec. 

00 

For the vortex "flow alone, it is found, that the kinetic 
energy/free stream velocity is 0.418 Newton sec. 

Hence, the energy change due to the acceleration of the 
particles is 0.53% of the total energy content of the vortex. 

It is felt that this is a small fraction of the vortex energy, 
and of itself, unlikely to strongly influence the vortex decay. 
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In view of this result, it is considered, that if water 
content in a vortex does indeed, influence vortex decay rates, 
then this must be enacted, through some mechanism other than 
energy absorption during particle acceleration. 
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VII. Discussion 


The study of particle motion within a vortex environment 
tended to indicate a general outward, migration of particles 
from the center of the vortex. Very heavy particles tend, to 
be more conscious of gravity and fall more directly through 
the vortex. Such particles are more rain-like and not the 
subject of the present study. 

The two-phase nature of the flow implies changes in the 
physical properties of the fluid.. These changes will manifest 
themselves as changes in the kinematic viscosity through both 
density and absolute viscosity variations. Methods of esti- 
mating these properties of two-phase mixtures are discussed 
in Two conventional expressions from^^are: 

|xt p = xp -w +(1_X ^g 

I - x , (1-x) 

P tp D w p g 


where x is the mass fraction of the liquid phase. Table 1 


shows the predicted variation of 


tp’ 


P t and v t with liquid 


mass fraction. The significant observation from this Table is 
that v^p decreases with an increase of water content. 

The influence of viscosity on the vortex decay was given 

(13) (17) 

by Squire , for turbulent flow and by McCormack et al v , 


for laminar flow. The latter result is: 

2 

= 1 + 


w 


( %ax ) 


4v ln2 
© 


v A 


( 11 ) 
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Using values from an example calculation of Ref. (17), equation 
(11) produces the vortex decay effects shown in Table 2. 

The calculations indicate that the molecular viscosity 
gives practically no vortex decay (in 20,000 ft) although 
what little trend there is gives a reduction in decay rate with 
increasing water content. The decay predicted, by Squire’s 
model for a turbulent vortex is shown in Figure 10. This 
figure demonstrates the effect of variations in eddy viscosity 
on the decay rate in an attempt to simulate the two-phase 
nature of the flow. Also shown in Figure (9) is experimental 

( 14 ■) 

data 'for a vortex with the same circulation as that adopted 
in the theoretical study. The large miss-match in slope may 
be due to the fact (as noted above for laminar vortices) that 
the theory is only valid very far downstream. 

Both theories predict large increases in vortex decay 
rate with increasing eddy viscosity, also, McCormack’s experimental 
data shows a significantly larger rate of decay than that calculated 
for laminar flow which is a result of the eddy viscosity of 
turbulence. Hence the presence of fog can still, conceivably, 
result in a more rapid, decay of the vortex if moisture serves 
to enlarge the eddy viscosity despite its effect on molecular 
viscosity. There is, however, very little information on the 
influence of water content on eddy viscosity and conflicting 
reports on two-phase mixing lengths give both larger and smaller 
values It is known that the liquid phase does significantly 
influence mean velocity profiles and this indicates that it may 
exert a strong effect on the eddy viscosity. How the droplets 

would interact with the turbulence eddy structure is not known. 

This problem remains to be resolved.. 
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The outward migration of droplets will have a tendency 
to cause a redistribution of the vortex velocity profile with 
the faster moving vortex core being retarded. It is possible 
that this could influence the mechanism of vortex bursting — 
particularly if the relationship between the axial and the 
asimuthal velocity components is distrubed.. Also the re- 
distribution of the droplets in the downstream direction results 
in a variation of the molecular viscosity across the vortex 
core. Further study is required, to understand, these effects on 
vortex bursting. Studies on vortex stability ' indicate 
that this phenomenon is closely related to the velocity gradients 
within the vortex. The above remarks strongly suggest that the 
water content in the vortex will go some way towards distorting 
the velocity profiles: and hence influence the flow stability. 

Another energy exchange mechanism may also be at work. 

The moisture saturated air in the higher pressure outer flow 
could suffer condensation and thus be warmed, while the low 
pressure core flow would be cooled, with attendent evaporation. 
This is an additional factor not included in the present model 
which may influence the vortex decay rate. 

It is believed that a considerably greater insight into 
the above mentioned effect of fog on trailing vortex decay 
could be obtained, from a more sophisticated, analysis of the 
Navier-Stokes equation than that reported, upon herein. 
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vxix. Concluding Comments 


The above exploratory remarks concerning the nature of 
vortex motion in an atmosphere containing fog have brought to 
light certain facets of the situation that would merit further 
deliberation. Among these we may mention the following: 

1) The results indicate that the kinetic energy transferred 
from the vortex to the particle during their acceleration 
is small. Thus water droplet acceleration is not likely to 
augment the rate of vortex decay. 

2) This conclusion does not exhaust the possibility that 
other mechanisms of vortex decay will not be accelerated 
by the presence of fog in the fluid medium. The two- 
phase character of the fog for example can conceivable 
enlarge the eddy viscosity, even though at the same time 
the molecular viscosity is decreased, and thus increase the 
decay rates or percipitate early vortex bursting. 

3) The method of the above analysis provides predictions 

of the velocity and particle trajectories subject to the 
assumptions of small particles and laminar inviscid vortex 
flow. It is suggested that the viscous term should be 
included and turbulence considered. 

4) The somewhat inconclusive findings of this tentative study 
are not presented to imply that fog will not enhance vortex 
decay, but rather, to illustrate that many unanswered, 
questions concerning this problem require further study 
for resolution. 
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Table 1 


Predicted Variation of Physical Properties 
With Mass Fraction of Water 


X 

P T p 

Hrpp 

7 TP 

0.01 

0. 70 

1.83 x 10~ 5 

2.63 x 10~ 5 

0.05 

3.19 

4.36 

1.37 

0. 10 

6.30 

7. 52 

1. 19 

0.20 

12. 50 

13.84 

1. 10 
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Table 2 


Calculated Vortex Decay rates using McCormack’s result (17) 



0.01 2.63 x 10~ 5 5.234 x 10" 8 -2.62 -5.23 -7.84 -10.46 


0.05 1.37 2.73 -1.36 -2.73 -4.09 - 5.45 

to 

^ 0.10 1.19 2.37 -1.19 -2.37 -3.56 - 4.74 
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Appendix Solution of the Vortex Equations 

It would, be most satisfactory if equation (4) of the 
main text were solved for a region of infinite radius. This 
solution can, of course, readily be written in terms of a 
Hankel transform, then if 

O 

A 1 


it follows that 


A(|x) = / r U (r) J (|ir) dr 
o 


U(B,z) = f |j.A((j,) J o (|j,r) exp (- \iz ) dp. A2 
o 


with similar expressions for v,f and. cp. 

If the form 

U 0 (r) = k 1 /(l+k 2 r 2 ) 

is taken for the initial profile, the equation A1 integrates 
to yield: k 


A(|-i) = 


1 K Q (n/k) 


with K q ( 5) the modified. Bessel function. 

Some algebra then gives an expression for U(r,z) in 
terms of transcendental functions. Specifically we have: 


2 2 - 


U(r,z)= k lV ^r 2 

I" 0 T 




with l = 2(l+i) 

Here r (5) and. F(n,m; -t; 5) denote the gamma and. hyper- 
geometric functions respectively. 

It can easily be demonstrated., however, that this series 
only converges for a modest range of r and z and is not 
suitable for a general study. 
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Attempts to solve the above integrals by means of Abel 
(15) 


transforms ' were not made. 

On the axis equation A2 integrates directly to give: 

5 . kz-1 


4 k x 


D(0 > z > - 3CTHS5 2 I •- -eir> 


A series expansion than takes the form 


4ki f o 

u(0 ’ z >- siiTkz ) 2 1 1 + ! 


kz-l 

kz+1 


•••) 


A3 
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Figure 2. Motion of Mean Radius Particle 
(Equation 2). 
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Figure 5. Distribution of Axial Velocity Component at Different Stations 
Along the Vortex. 





WORK ft 




34 





